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^3 , Abstract. Information about the X-ray luminosity of the supermassive black hole located at the Galactic center 

(GC), Sgr A*, and its temporal variations in the past is imprinted in the scattered emission observed today in 
the direction towards giant molecular clouds (GMCs) located in our Galaxy. Due to light travel time effects these 
clouds probe the activity of Sgr A* at different times in the past depending on their position relative to the 
GC and the observer. In this paper we combine results of recent ASCA observations along the Galactic plane, 
, providing upper limits for the scattered flux in the 4-10 keV range produced in a given direction, with data from 

^ ■ CO surveys of the same regions. These CO surveys map the position and mass of the molecular gas which the 

' GMCs are made up of. Demanding the scattered flux to be not larger than the observed one, this data enables 

^\ , us to derive upper limits for the 4-10 keV luminosity of Sgr A* at certain times during the last 40 000 years 

' down to about 8 x 10'"'ergs~^. At other times the limits are less tight, of the order of 10*^ — lO^^ergs"^. For 

two periods of time of about 2000 and 4000 years duration 8000 and 14000 years ago the currently available 
^ ■ CO data is insensitive to any enhanced activity of the GC. Flares lasting longer than 3000 years flU these time 

gaps and therefore can be excluded to have occured during the last 40 000 years with a luminosity larger than a 
few lO^^ergs-V The more extended and continuous HI distribution in the Galactic disk, which also scatters the 

■ radiation emitted by Sgr A*, allows us to extend the time coverage further into the past, back to about 110 000 
years, albeit the limits are becoming less tight. We thereby can rule out a long term X-ray activity phase of Sgr 
A* at one per cent of its Eddington level ending less than about 80 000 years ago. The limits presented in this 

5^ ' paper can be improved by observations of emission in the fluorescent iron Ka-line. We study the feasibility of 

^ ' these methods to investigate past nuclear activity in other spiral galaxies observed with the angular resolution of 

. X-ray telescopes like Chandra and XMM-Newton. 
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1. Introduction from CO observations. Sunyaev at al. (1993) proposed that 

the diffuse, hard X-ray emission observed today in the di- 

High spatial resolution observations of the stars in the ^^^^j^^ g^^^^ molecular clouds (GMCs) in the central 

vicinity of the Galactic center (GC) strongly suggest the ^^1^^^ especially Sgr B2 is radiation 

presence of a supermassive black hole of mass 2.6 x 10^ Mq ^^^^^^^ scattered into our hue of 

(Genzel et al. 1997; Chez et al. 1998). Recent Chandra ^^^^^ molecular gas. They predicted the presence of 

measurements (Baganoff et al. 2001) detect the associated ^ ^^^^^^ -^^^ fluorescent K„-line at 6.4 keV with an equiv- 

X-ray source Sgr A* at a mere 2 x 10^^ erg s , more than ^^^^^ ^j^^j^ ^j^^ ^^^^^ ^ r^^-^ ^ signature of 

ten orders of magnitude below the Eddington luminosity g.^ttering of hard X-ray photons on neutral matter when 

of 3 X 10^4 erg s ^ for a black hole of this mass. This rela- ^^-^^^^ ^^^^^^ ^^^^^^ ^^^^ contribute to the ob- 

tive dimness being remarkable in itself raises the question ^^^^^^ emission, by either falling outside of the observing 

whether Sgr A* was equally dim in the past. ^^^^ ^^-^^ currently not active (Vainshtein & Sunyaev 

Granat spacecraft observations of the region surround- iggQ) 
ing the GC have revealed a close resemblance between the 

morphology of the diffuse X-ray emission above 10 keV o^der for the mass of molecular gas (~ 6 x 10^ Mq) 

and the spatial distribution of molecular gas as inferred present in Sgr B2 to produce the observed flux the X- 

ray luminosity of Sgr A* must have been of the order of 

Send offprint requests to: conrad@mpa-garching.mpg.de 2 X lO^^ergs^^ about 400 years ago. This time delay corre- 
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spends to the different path lengths for photons received 
directly or taking the detour via Sgr B2, which lies at a 
projected distance of ~ 120 pc from the GC. 

Evidence for the scattering scenario has been discov- 
ered with the spectroscopic imaging capabilities of the 
ASCA satellite (Koyama et al. 1996; Sunyaev et al. 1998; 
Mnrakarni et al. 2000). The X-ray spectrum of Sgr B2 ex- 
hibits strong emission in the iron Ko,-line at 6.4 keV con- 
firming the prediction of Sunyaev et al. (1993). Murakami 
et al. (2001A) have recently observed another GMC close 
to the GC, Sgr C, to be bright in the K^-line, strength- 
ening the case for the scattering scenario. Theses authors 
call Sgr B2 and Sgr C "X-ray reflection nebulae" . Recent 
Chandra observations (Murakami et al. 2001B) confirm 
the ASCA results about Sgr B2 and the nature of the 6.4 
keV Kc-line emission. 

The application of the method successfully used in the 
above described cases to constrain the luminosity of Sgr 
A* in the past is not restricted to GMCs close to the GC. 
For it to provide nontrivial constraints upon the strength 
of the activity of the GC one needs a scattering mate- 
rial with a sufficient optical depth and an estimate of the 
scattered flux produced by this material. In the case of 
our Galaxy two obvious choices which fulfill these criteria 
are interstellar molecular and atomic hydrogen: 

a) CO surveys covering the Galactic plane have shown 
that about 10^ Mq of H2 reside in GMCs on orbits 
between the GC and the sun with a peak mass sur- 
face density around a galacto centric radius between 4 
and 6 kpc, the so called molecular ring (e. g. Dame 
1993). The time delays, the difference of the arrival 
times between the primary and the scattered radia- 
tion, for clouds located in the molecular ring range up 
to about 40 000 years. 

In this paper we will compute the scattered flux pro- 
duced by a GMC in response to a strong flare of Sgr 
A* . This flux is detectable by current X-ray telescopes 
even if it was relatively short, with a duration of the 
order of years. Turning this argument around one in- 
fers from the observed relative dimness of GMCs today 
that there could not have been a strong flare of Sgr A* 
at those times in the past, which correspond to the 
time delays of massive GMCs in our Galaxy. 
The good temporal "resolution" with an uncertainty 
of the occurrence of a flare of the order of years comes 
about because the mass of a GMC is rather concen- 
trated with a very dense core. We are therefore in prin- 
ciple able to pinpoint past flares of Sgr A* during the 
last 40 000 years very well if there exists a GMC with 
a corresponding time delay. 

The drawback of this "accurate" timing is that because 
there only exist a finite number of GMCs in our Galaxy 
it is not possible to cover the whole time interval from 
40 000 years ago, the time delay of the most distant 
clouds, until "today". For the analysis presented in 
this paper it was only possible to use the CO data 
from the first Galactic quadrant. This means that we 



are missing all the GMCs, which are located in the 
fourth Galactic quadrant and which should be equally 
numerous. This additional data would probably im- 
prove the "filling factor" of our sample but still leave 
several time gaps and not extend the time coverage 
beyond 40 000 years. The possibility thus always ex- 
ists that a short flare of Sgr A* occurcd during the last 
40 000 years unnoticed by us, because there is no GMC 
which could respond to this flare. If we take only the 
available CO data into account the coverage becomes 
complete for flares longer than about 3000 years and 
we can thereby rule out any flare brighter than a few 
10^^ erg s^^ during the last 40 000 years of this or a 
longer duration, 
b) The HI phase of the interstellar medium (ISM) is more 
diffusely distributed than H2. It extends out to about 
16 kpc with an approximately constant mass surface 
density of about 4M0pc~^, corresponding to a total 
mass of the order of 3 x 10^ Mq (e. g. Dame 1993). 
This implies that the total mass of neutral hydrogen 
in a field of view of the ASCA instrument is larger than 
the mass of any GMC even those near the GC. For a 
Galactic longitude of 10° it is about 2 x 10^ M©. 
Due to its more continuous distribution the response 
of the HI distribution to flares of Sgr A* will be spread 
out over a long time interval. Therefore we are unable 
to locate the time when a flare happened. A strong 
flare longer ago can produce the same response as a 
weak flare not so long ago. Nevertheless due to the 
larger extent of HI in the Galactic disk we can ex- 
tend the time coverage back to about 110 000 years. 
We will show that for example the switch off case from 
1% of the Eddington luminosity of Sgr A* (the whole 
Galactic disk is illuminated by Sgr A* , which suddenly 
"switches off" at a time to) can be ruled out to have 
ended less than about 80 000 years ago. 

The ASCA Galactic Plane Survey (Sugizaki et 
al. 2001), provides upper limits for the diffuse X-ray emis- 
sion produced by scattering along certain directions. In 
this paper we compute the expected scattered flux for an 
assumed Eddington luminosity of 3 x lO^^ergs"^ in the 
4-10 keV range. Although this value is rather large the 
results can be rescaled to any input luminosity because 
the problem is linear. It is obvious that the ratio of the 
computed scattered flux to the upper limit for the scat- 
tered flux given by ASCA observations, / = i^sc/-Pobsj 
gives us an estimate of the average luminosity of the GC, 
^GC ^ .f "^^Edd, at a time which corresponds to the time 
delay of the observed GMC. 

This paper is structured as follows. In the following 
section wc recall the basic ideas and relations of X-ray 
archaeology. We describe how we extracted the CO and 
X-ray data necessary for our purposes from the literature. 
In Sec. 3 we will derive constraints for the past X-ray lumi- 
nosity of Sgr A* using a sample of GMCs in the Galaxy. In 
Sec. 4 we will describe the limits the Galactic HI distribu- 
tion can provide. The possibility to use the same methods 
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Fig. 1. A sketch of the scenario studied projected onto the 
Galactic plane. A time to ago the luminosity of Sgr A* 
(marked by a cross) flared. This enhanced activity lasted 
for a period A. Relative to the earth bound observer (small 
circle), observing at a time t and orbiting Sgr A* on the 
dotted circle, the GMCs that lie between the two ellipsoids 
corresponding to the beginning and the end of the flare are 
exposed to this enhanced emission. For other GMCs the 
period of enhanced irradiation already has passed or has 
not started yet. The boundaries of an observational field of 
view are sketched by dashed lines. Inside this field of view 
one cloud has already faded and another one is partially 
covered by the flare. Two possible photon paths between 
Sgr A* and the observer are marked by dot-dashed lines 
(one of them covered halfway by the z-axis). The ellipsoids 
are rotationally symmetric about the z-axis. 



to study past nuclear activity in other spiral galaxies is 
studied in Sec. 5. In Sec. 6 we close with a discussion of 
our results. 



2. Method 

2.1. X-ray Archaeology 

The scenario we have in mind is sketched in Fig. 1. An 
outburst of duration A has occured at the GC a time to 
ago, as measured by the earth bound observer. As time 
evolves this signal moves radially outwards from the GC 



and starts illuminating GMCs at larger and larger galac- 
tocentric distances. These clouds scatter a fraction of the 
incident radiation into the direction of the observer. 

The loci of scattering sites giving a fixed time delay 
form an ellipsoid with foci at the source and the observer. 
Since the system is symmetric around the line connecting 
the GC and the observer this ellipsoid is fully described 
by two parameters which depend upon time past since 
the detection of the burst t — to, where t is the time of 
observation. In this section we will assume from now on 
to = to simplify the following formulae. 

Choosing a cartesian coordinate system (see Fig. 1) 
with the GC and the Sun at z — Ro/2 and z = —Ro/2 
respectively, the ellipsoid can be expressed as 



y 



62 



1, 



(1) 



where the lengths of the major and minor axis are a = 
{ct + Ro)/2 and b = y/ict)'^ + 2ctRo/2, with c the speed 
of light and Ro the distance to the GC. 

An alternative description of the ellipsoid can be given 
in polar coordinates R, (j) (defined in Fig. 1): 



1 



e cos ( 



with p = b^/a and e = Ro/{2a) 
of the observer at infinity {Rq — 
ellipsoid becomes a paraboloid 



i?(0) 



ct 



1 -I- cos ( 



(2) 

Note that in the limit 
oo,p ^ ct,e 1) this 

(3) 



This approximation is valid for GMCs at the GC {Ro ^ 
ct) and has been studied in the context of Sgr B2 by 
Sunyaev & Churazov (1998).0 In cartesian coordinates the 
paraboloid is described by: 



1 / X- 
z — — \ ct 
2 



ct 



Ro 
'Y' 



(4) 



The velocities of the surfaces of constant time delay 
along the major and minor axes of the ellipsoid are ob- 
tained by differentiation of the formulae for a and b. The 
major axis increases at a constant speed of c/2. The ve- 
locity along the minor axes depends upon time as follows 



ct + Ro 



2 y^{ct)^ + 2ctRo 



(5) 



At times small compared to Ro/c this velocity is larger 
than the velocity of light and approaches the value c/2 

^ The study of light echoes has a long history, starting with 
an explanation of observations of Nova Persei 1901 by Couderc 
(1939). Since then it has been applied to several astrophys- 
ical objects, e. g. supernova lightcurves (Morrison & Sartori 
1969), scattering by interstellar dust (Alcock & Hatchett 1978), 
scattering of an isotropic (Sunyaev 1982) or anisotropic radio 
source (Gilfanov et al. 1987) by the thermal gas in a clus- 
ter of galaxies, reverberation mapping in AGN (e. g. Peterson 
1993) and Compton echoes from gamma-ray bursts (Madau et 
al. 2000) to name a few. 
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for ct ^ Ro- Inside the Galaxy, where most of the GMCs 
lie inside the solar circle {R < Rq), this velocity is always 
larger than (2/\/3)(c/2) « 0.58c. 

The velocity in a general dir(X'tion (f> is given by the 
temporal derivative of Eq. 2 which yields 



c c^V + 2i?o(ct + i?o)(l + cos( 
2 (ct + i?o(l + cos0))2 



(6) 



In the paraboloid approximation this becomes, by either 
differentiating Eq. 3 with respect to t or by taking the 
limit Ro > ct of Eq. 6, 



R{<t>) 



1 + cos ( 



(7) 



The apparent velocity is given by the projection onto the 
plane of the sky, which reads 



csm( 



^app 



1 + cos ( 



(8) 



This equation is equivalent to the formula for apparent 
superluminal motion of matter having velocities close to 
the speed of light, which is normally expressed in terms of 
the angle 9 = tt — (f>. The formula for the apparent velocity 
in the Galactic plane in the ellipsoid case is given by 



^app 



c cH'^ + 2Ro{ct + i?o)(l + cos( 

2 (ct + i?o(l + COS0))2 



sm{(P-l). (9) 



One obtains an additional dependence upon I, because 
the projection angle for a given direction (f> depends upon 
galactocentric distance R and therefore upon time. 

The average time delay of photons scattered by a cloud 
with galactocentric radius R and distance from the ob- 
server d is given by: 



6={R + d- Ro)/c. 



(10) 



It is obvious that for clouds with a constant galactocentric 
distance R those, which lie behind the GC, produce the 
largest time delay 5m^^{R = const.) = 2R/c. For clouds 
with a distance of 6 kpc to the GC, the edge of the molecu- 
lar ring, the largest time delay is about 40 000 years. The 
more diffuse HI phase of the interstellar medium (ISM) 
extends up to radii of about 16 kpc (Dame 1993). Matter 
at these distances can in principle probe the activity of 
the GC about 110 000 years ago. 

2.2. GMCs in the Galactic disk 

The results of several CO surveys have been published 
over the last decades. They map the CO brightness tem- 
perature as a function of Galactic longitude I, Galactic 
latitude b and radial velocity component Vr, differing in 
resolution and area surveyed. 

Solomon et al. (1987, SRBY87) have published a list 
(see Table 1 in their paper) of 273 molecular clouds de- 
tected with the Massachusetts-Stony Brook CO Galactic 
Plane Survey of the first Galactic quadrant. These clouds 
were extracted from the three dimensional {l,b,Vr) data 
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Fig. 2. The spatial distribution of GMCs selected from 
SRBY87 for this work. The GMCs cover a range of 
Galactic longitudes of 8° < / < 60°. The position of 
the GC is marked by a cross. The position of the sun 
(marked by a small circle) is at coordinates (—8.5,0.0). 
For Galactic longitudes c;overed by ASCA observations 
the solar circle is marked by solid arcs. For Galactic longi- 
tudes falling outside of the ASCA survey the solar cir- 
cle is marked by a dotted line. The different symbols 
denote different cloud masses: M^i < 5x10^ Mq (di- 
amonds), 5 X 10^ M© < Mci < 10^ Mq (triangles) and 
Mci > 10^ Mq (squares). 



cube by locating surfaces of equal brightness tempera- 
ture. The distances were computed from the radial ve- 
locity component of the centroids under the usual as- 
sumption (see e.g. Binney & Merrifield 1998) that clouds 
move on circular orbits with a constant circular velocity 
of Oo = 220 km s"^ 

The assumption of circular orbits with a constant ro- 
tational velocity is a fair approximation for clouds inside 
of the solar circle but breaks down at small galactocentric 
radii because the orbits are probably highly non circular 
in this region. For GMCs inside the solar circle a ambigu- 
ity is inherent to this method of distance determination. 
A cloud with a given measured radial velocity component 
can either lie on the near or the far side of the circular orbit 
it is assumed to travel around the GC. SRBY87 were able 
to overcome this ambiguity by several methods described 
in their paper. 

This data set is not the most recent one (for more 
recent surveys of this area see, e.g. Dame et al. (2001) 
and Lee et al. (2001)) but to our knowledge up to today 
the only one where mass estimates of a large number of 
individual GMCs have been presented. 

Not all of the clouds listed in SRBY87 fall into one 
of the fields of view of the ASCA Galactic Plane Survey. 
These are clouds located at longitudes larger than 60° or 
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^ Time delay 5 and time for the ellipsoid to cross the cloud td in units of years 

Mass of cloud in units of 10* M© 
" Scattered flux in the 4-10 kcV range in units of crgcm"^ s"! for the whole cloud covered by a flare 

Flux in the 4-10 keV range observed by ASCA in units of ergcm""^ s^i 
" Upper limit for the luminosity of Sgr A* in the 4-10 keV range in units of ergs-i a time S ago 

Table 1. Selected parameters (number of cloud in SRBY87 sample, Galactic longitude, Galactic latitude, galactocentric 
distance, distance, time delay, mass, radius, time for ellipsoid to cross the cloud, scattered flux, observed flux and upper 
limit for the luminosity of Sgr A*) of twenty GMCs from our sample. 



at latitudes larger than 0.3°. Combining the two data sets 
leaves one with 96 clouds out of 273. The spatial distri- 
bution of these clouds in the Galactic plane is plotted in 
Fig. 2, symbolised according to their mass (for the dis- 
tribution of the full set of clouds see Solomon & Rivolo 
1989). Although spiral structure might be marginally dis- 
cernible in this picture, the concentration of GMCs along 
the molecular ring is clearly visible. 

The cloud properties needed for this work are the 
Galactic longitude the galactocentric radius R, the dis- 
tance from the observer d, the radius and the mass 
Mci of a cloud, which was computed with the observed 
velocity dispersion and should thus not be affected by re- 
calibrations of the I(C0)/N(H2) ratio. For the subsample 
of clouds we are working with these parameters fall in the 
ranges 1.9 x lO^M© < < 3.8 x lO^ Mq, 2.5 < < 47 
pc and 3 x lO^years < ^ < 4 x lO^years. Note that in 
SRBY87 a distance of the sun to the GC of 10 kpc was 
assumed. We have rescaled the distances and derived cloud 
properties to Rq = 8.5 kpc. In Table 1 (columns 1-8) we 
list the important parameters for twenty clouds out of 
our sample. Note that the GMCs in SRBY87 have been 
ordered and numbered according to their Galactic longi- 
tude. 



2.3. Scattered flux 

The ASCA observations of the Galactic plane (6 w 0°) in 
the longitude range —45° < / < 63° have been compiled 
and analysed by Sugizaki et al. (2001). In Fig. 2 we have 
tried to illustrate the range of longitudes covered by the 
ASCA Galactic Plane Survey by plotting solid arcs at a 
distance -Ro from the GC, where observations have been 
made;. It is obvious that we are missing the CO data from 
the fourth Galactic quadrant. As mentioned before sur- 
veys of this area exist (e. g. Dame et al. 2001) but to our 
knowledge this data has not yet been fully analysed. The 
distribution of GMCs should be symmetric between the 
Galactic quadrants I and IV. We should therefore expect 
a similar number of GMCs in the fourth Galactic quad- 
rant. 

The ASCA observations have been presented in three 
energy bands. Most interesting for this work because least 
affected by photoelectric absorption is the X-ray surface 
brightness in the hardest, the 4-10 keV, spectral band. 
The corresponding fluxes at all longitudes are all of a sim- 
ilar order of magnitude. Not even one cloud seems to be 
especially bright "at the moment" . 

Even if the diffuse X-ray flux observed in one field 
of view originates from a different physical process than 
scattering, it is still an upper limit for the flux produced 
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through scattering by gas distributed along this beam: 

-^sc — -^obs* 

The clouds in our sample have optical depths for 
Thomson scattering up to tt ~ lO^""^. Wo therefore c;on- 
sider scattering in the optically thin limit. Our treatment 
of absorption is described in Sec. 2. 3. 1. Nearly the whole 
Galaxy is transparent to X-rays with energies above 4 keV, 
except some regions around the GC. But even here op- 
tically thick GMCs do not shield the whole disk of the 
Galaxy from the illumination by Sgr A*. The covering 
fraction is probably rather small. 

In the optically thin case the scattered surface bright- 
ness in a given direction at energy E and time of observa- 
tion t is given by 



Isc{E,t) 



At) 



n{s) 



Lgc{E, ^rct) 

47ri?(s)2 



da 



ds, (11) 



„i„(t) 

where Smin and Smax are the boundaries of integration, 
n(s) is the number density of the scattering gas at a dis- 
tance s from the observer, tret =t— {R{s) + s)/c the "re- 
tarded" time, R the distance to the GC and (da/dil) the 
differential scattering cross-section. The time dependent 
integration limits are the points of intersection of the line 
of sight with the ellipsoids corresponding to the beginning 
and the end of a flare. 

The scattered flux is given by the integration of the 
surface brightness over the field of view of the ASCA in- 
strument. We take into account the cases when parts of a 
cloud fall outside the field of view, as sketched in Fig. 1, 
or when the field of view is larger than a cloud: 



/ 



F,^{E,t) = / I,c{E,t)dQ. 



(12) 



We neglect here the iron K„-linc production and the cor- 
responding iron K-edge absorption. The detailed ASCA 
spectral data are not yet publicly available. XMM-Newton 
will provide much improved data including K„-linc limits, 
which will improve the limits for the X-ray luminosity of 
Sgr A* derived in the following sections. 

In the optical thin case the scattered luminosity is pro- 
portional to the mass of the illuminated scattering mate- 
rial. An order of magnitude estimate of the above integrals 
to demonstrate the scaling yields 

This formula neglects the spatial extension of the scat- 
tering gas and the angular dependence of the differential 
scattering cross-section and assumes that the whole cloud 
is covered by the flare, implying that all of the mass con- 
tributes to the scattered emission. For typical cloud pa- 
rameters one obtains a scattered luminosity of 



L 



Edd 



10« Mr.- 



R 



which for Lgc ~ -^Edd is comparable to the luminosity of 
the brightest X-ray binaries in the Galaxy, thereby cer- 
tainly contradicting the observations. 



If the flare duration A is shorter than the time for the 
ellipsoid to sweep across the cloud one has to multiply 
Eq. 14 by an additional factor 



(15) 



where fdi is the velocity of the ellipsoid across the cloud, 
which can be superluminal depending on the relative ori- 
entation of the cloud with respect to the GC and the ob- 
server. This additional factor takes into account that only 
a fraction ~ A/(2rci/i'en) of the mass of the cloud is illu- 
minated. 



2.3.1. Absorption 

As X-ray photons of energy E traverse the ISM a frac- 
tion Tph{E) of them is absorbed by the heavy elements 
of the interstellar gas, again assuming optically thin con- 
ditions. Therefore a GMC will be exposed to a flux 
(Lgc(£, tret)/47ri?2) exp[-r^h(^)]> where t'^^{E) is the op- 
tical depth for photoelectric absorption at energy E be- 
tween the GC and the GMC. The observer will measure 
a flux Fsc(£^, i) exp[— Tph(£')]. The optical depth due to 
photoelectric absorption is proportional to the hydrogen 
column density A^h 



rph(£;) = N^<Jp^{E), 



(16) 



where Cph {E) is the photoelectric absorption cross-section 
per hydrogen atom, which we compute with an analytic 
fit by Morrison & McCammon (1983). 

We estimate A^h with the average mass surface density 
of a GMC given by Sci ~ -^ci/^'ci''') which leads to A^h ~ 
(3/4)(Sci/mp). Since we have no better knowledge about 
the hydrogen column density between the GC and a cloud, 
we furthermore assume = N^. 

2.3.2. Differential scattering cross-section for a 
mixture of H2 and He 

In the following sections we will assume that the GMCs 
in our sample are composed of 70% molecular hydrogen 
and 30% helium. It is known (e. g. Sunyaev & Churazov 
1996) that in the energy range considered scattering by 
the bound electrons of molecular hydrogen and helium 
is important, because for small scattering angles the dif- 
ferential scattering cross-section per electron for Raylcigh 
scattering on molecular hydrogen and helium can become 
two times larger than the one for hydrogen atoms, which 
at these energies is identical to the Thomson differential 
cross-section. The angular dependence for Thomson scat- 
tering is proportional to (1 -|- cos^^sc); where is the 
scattering angle of the photon, which is approximately 
given by the angle a defined in Fig. 1. 

In Fig. 3 we plot the scattered flux for a single cloud 
of 10® Mq composed of pure molecular hydrogen, a mix- 
ture of 70% molecular hydrogen and 30% helium and pure 
atomic hydrogen as a function of the angle a {a = w cor- 
responds to the case of forward scattering). We chose a 
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Fig. 3. The scattered flux in the 4-10 keV as a function 

of the angle a (defined in Fig. 1) for a cloud of lO*' Mq at 
5 kpc distance from the GC and the observer fully illumi- 
nated by an Eddington flare of Sgr A*. Note that a = n 
corresponds to the case of forward scattering. The lines 
mark different mass compositions of the cloud: 70% H2 + 
30% He (solid line), atomic hydrogen (dashed) and molec- 
ular hydrogen (dotted). 

galactocentric distance and a distance to the observer of 
5 kpc. We used the diff'erential scattering cross-section for 
scattering by atomic, molecular hydrogen and helium as 
a function of photon energy and scattering angle as pre- 
sented in Sunyaev et al. (1999). 

The scattered flux for the compositions including 
molecular hydrogen are larger than the one for atomic 
hydrogen because for small scattering angles (a < tt) the 
two electrons present in molcciilar hydrogen and helium 
scatter coherently. The scattered flux for the chosen mix- 
ture of molecular hydrogen and helium is smaller than the 
one for pure molecular hydrogen because in this case there 
are only .85 electrons per nucleon compared to 1 electron 
per nucleon for pure hydrogen, atomic or molecular. 

3. Limits on the X-ray luminosity of Sgr A'^ from 
Giant Molecular Clouds 

There exist several characteristic variability timescales for 
a supermassive black hole such as Sgr A*. The short- 
est possible timcscalc is related to the linear size of the 
region, where the bulk of energy is released in the ac- 
cretion disk around a black hole (Shakura & Sunyaev 
1973): fbh lORg/c - 300 s for Sgr A*, where i?g is the 
Schwarzschild radius of the black hole. Tidal disruption 
events are believed to produce an Eddington level flare 
over a period ttidai ^ years (Rccs 1988). The lifetimes of 
active galactic nuclei (AGN) are believed to be of the or- 
der of 10^ — 10^ years. The detailed temporal evolution 
of the luminosity of Sgr A* is probably rather complex, 
possibly varying on all of the above timescales. 

In this study it is our aim to constrain the past lumi- 
nosity of Sgr A* from information about the reprocessing 



system, the spatially distributed GMCs in the Galactic 
disk, and the output signal at one point in time, the 
scattered flux observed today. This is reminiscent of the 
methods used in reverberation mapping in AGN (see c. g. 
Peterson 1993). There one tries to gain information about 
the geometry and kinematics of the emission-line region 
(the reprocessing system) by monitoring the continuum 
variations (the input) and the resulting emission-line re- 
sponse (the output). Mathematically this is described by 
a transfer function, which applied to our problem can be 
written as: 

/oo 
■^{t')LGc{t-t')dt' (17) 
-00 

The scattered flux observed at a time t is the response of 
our system to activity of the GC in the past. In reverbera- 
tion mapping in AGN one determines the transfer function 
using Fourier transforms of the input and the output over 
a certain period of time. In our case this is not possible 
because the scattered flux or its upper limit is known only 
at one moment in time. 

If the cloud had no spatial extension, a "point mass" , 
the transfer function would be just a (5-function with an 
amplitude proportional to the mass of the cloud. The 
transfer function for several clouds with different time de- 
lays would be a sum of (5-functions. In these cases it would 
be possible to ascribe an unique X-ray luminosity of Sgr 
A* to an observed scattered flux. 

Due to the spatial extension of a GMC with an internal 
density distribution this ^-function is smeared out in time. 
This means that the scattered flux observed at one mo- 
ment in time depends upon the luminosity of the GC over 
a certain period of time in the past. The correspondence 
between the observed scattered flux and the luminosity of 
the GC is therefore not one to one anymore. 

What one can do in this case to constrain the luminos- 
ity of Sgr A* is to assume some input, a "test function", 
and compare the resulting output with the observations. 
Should the output be inconsistent with the observed flux 
then the initial input must be wrong. 

We assume the temporal and spectral behaviour of the 
luminosity of Sgr A* to be of the following form 

LGc(£;,i) = |o°^^°/^^ else"*"*"" ^^^^ 

with ton = to — A/2 and tos = to + A/2. The choice of the 
spectral behaviour implies a photon index of 2, which is 
representative of observed AGN spectra in the 4-10 keV 
range. Note that here we define for the sake of symmetry to 
as the time half way through the flare, whereas in Sec. 2. 1. 
to was deflned as the beginning of the flare. 

Different inputs can produce the same output because 
there is as mentioned above no one to one relation between 
the luminosity and the scattered flux. It is therefore im- 
possible to constrain the past luminosity at one moment in 
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time. The important quantity is thus actually the emitted 
energy over a period T 



Q 



LGc{t')dt', 



(19) 



where T is given by the minimum of the flare duration 
and the time for the ellipsoid to cross the cloud. The lim- 
its we are deriving in the following sections are therefore 
constraints upon a time averaged luminosity of Sgr A*, 
viz. Lgc < imax with Lqc = Q/T. 

3.1. Response of a single cloud 

Observations of molecular tracers of H2 (see e. g. Blitz & 
Williams 1999) show that GMCs are far from being spher- 
ically symmetric objects. They are highly structured being 
made out of clumps with high density variations from one 
point to another. Due to their amorphous nature without 
a clear central peak the density profile is often difficult 
to define observationally. Another difficulty is that differ- 
ent tracers of H2 turn optically thick at different densities. 
Nevertheless Williams et. al (1995) were able to determine 
a density distribution of the form p cx for dense clumps 
in the Rosette molecular cloud. We therefore adopt this 
radial dependence as our profile A: 



Pci(r-) 







r > Tel. 



(20) 



The normalisation po is chosen in such a way, that the 
mass contained inside the radius of the cloud corresponds 
to the one given in SRBY87. The cutoff is required to 
obtain a finite mass. 

Besides the response of the cloud with density profile A 
we have also computed the response for a different profile 
to study the influence of the mass distribution inside a 
cloud upon our results. We choose as our profile B the 
following one: 



pair) 



Po 



(1 + (r/ro)2)2 



(21) 



This profile has a constant core density up to a radius of 
the order Tq, which we choose to be O.lrd. Most of the 
mass of this profile is provided by matter located at these 
radii. 

Having specified the input (Eq. 18) and an internal 
density distribution for a cloud (Eqs. 20, 21) we can inte- 
grate Eqs. 11 and 12 numerically to compute the scattered 
flux produced by a cloud out of our sample. We chose cloud 
116 from the SRBY87 sample (see Table 1) and computed 
its response to Eddington flares of Sgr A* of durations of 
1, 10, 100 and 300 years as a function of time. Fig. 4 dis- 
plays the results of these computations .0 

^ Note that Fig. 4 (and the following figures depicting a de- 
pendence upon t — to) can be read in two ways: 

1. Regarding t, the time of observation, as the independent 
variable and keeping to, the time half way through the flare, 
fixed (e. g. to = ~S) the different curves mark the scattered 




Fig. 4. The scattered flux observable at a time t produced 
by a single cloud (number 116 in Table 1) in response 
to an Eddington flare of Sgr A* of 1, 10, 100 and 300 
years duration (from bottom to top) occuring around a 
time to for two different density profiles of the cloud given 
by Eq. 20 (solid hues) and Eq. 21 (dashed lines). The 
observational limit for the scattered flux produced in the 
field of view the cloud is located in is marked by a dotted 
line. The time delay S for this cloud is 20 764 years, the 
radius of the cloud is 13.4 pc and the time for the ellipsoid 
to cross the cloud about 126 years, which corresponds to 
a cloud crossing velocity of Weii ~ 0.8c (for the other cloud 
parameters see Table 1). 



It shows that an Eddington flare of the GC lasting 
for one year 20 764 years ago, the time delay S of cloud 
116, would produce "today" a scattered flux of about 2 x 
10"^" erg cm~2 s~^ for profile A and 4x lO"""^^ ergcm"^ s~^ 
for profile B in the ASCA field of view cloud 116 is located 
in. Since the observational limit for this field of view is 
below these values, thus falsifying the initial assumption, 
there could not have been an Eddington flare of Sgr A* of 
a duration of one year or longer 20 764 years ago. A one 
year flare occuring 100 years earlier or later can not be 
excluded by cloud 116, because the scattered fluxes would 
be below the observational limits. 

Because especially for short flares the shape of the 
lightcurve mainly depends upon how the mass is dis- 
tributed inside a cloud the responses for the density pro- 
files A and B differ. 



flux as a function of observing time t in response to an fiare 
of Sgr A* at to — ~S. One tracks the response of a single 
fiare as time evolves. 

Treating to as the independent variable and keeping t, the 
time of observation, fixed (e. g. t = 0) the different curves 
mark the scattered fiux observable at a time f = as a 
function of the time half way through the flare to. At one 
moment in time one observes the response of fiares occuring 
at different times to. This case applies to what is done in 
this paper to constrain the past luminosity of Sgr A*. 
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The response of a cloiid with a density profile A to the 
one year flare illustrates that a cloud whose mass is cen- 
trally concentrated will produce a light curve with a very 
sharp maximum. This comes about because for short flare 
durations the space between the two ellipsoids covers only 
a small part of the cloud. This illuminated region sweeps 
along the cloud and produces a maximum flux when it cov- 
ers the most dense, central part of the cloud. Increasing 
the flare duration one encloses more and more mass be- 
tween the two ellipsoids, corresponding to the beginning 
and the end of the flare. This increases the scattered flux, 
which is proportional to the mass, and averages out den- 
sity contrasts at different radii. The central spike produced 
by the one year flare therefore vanishes for longer flare du- 
rations. The steep decline of the scattered flux observable 
for density profile A comes about because of the cutoff of 
the profile at r = rd- 

For density profile B, which has a constant density 
inside Tq, the main contribution to the total mass of the 
cloud is provided by matter at radii r ~ Tq. Because the 
one year flare covers a region smaller than this size, there 
is more mass probed by the one year flare in the central 
region of profile A than profile B. The largest scattered 
flux produced for a one year flare by a cloud with density 
profile A is therefore larger than the one for a cloud with 
density profile B. For fiares of longer duration larger parts 
of the cloud will be covered and the fluxes will reach a 
similar level, because both density profiles are normalised 
in such a way that the total mass contained inside of a 
radius rd is equal for both profiles. 

If the duration of the flare increases, the whole cloud 
will be covered at one point. From this point on an increase 
of flare duration will no more produce a larger flux, be- 
cause all of the cloud is already enclosed between the two 
ellipsoids. This "saturation" time tci, the flare duration 
for which the maximum flux is achieved or equivalently 
the time for the ellipsoid to sweep across the cloud, de- 
pends obviously on the radius of the cloud but also on the 
position of the cloud relative to the GC and the observer, 
which determines time travel effects and the different ve- 
locities of the ellipsoid sweeping across the cloud. For ex- 
ample for the extreme case, when the cloud lies exactly 
behind the source the time it takes an ellipsoid to cross 
the cloud is 4rci/c. On the contrary for the case when 
the cloud lies half way between the GC and the observer 
at a small Galactic longitude, the time to sweep across 
the cloud is of the order of {rd/ Ro)rci/c, which is obvi- 
ously much less than the value above. Thus if the cloud 
is located closer to the observer than the source of the 
flare, the whole cloud will be scanned much quicker than 
the light crossing time of the cloud. For the cloud behind 
the source the temporal evolution will be slowest. For the 
GMCs in our sample the values lie somewhere between 
these two extreme cases. 

We have computed lightcurves similar to Fig. 4 for 
all the GMCs in our sample. Therewith we are able to 
estimate the saturation time and the maximum scattered 
flux for them. The saturation time td of a cloud and its 
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radius yield an estimate for the velocity of the ellipsoid 
sweeping across the CMC: Vdi ~ {2rd/c)/td- In Table 1, 
column 9 we have listed the saturation times for twenty 
GMCs from our sample. The cloud crossing velocities for 
these clouds lie in the range 0.6c < Vdi ^ 3.2c. 

The maximum scattered flux produced by a cloud 
yields an upper limit for the luminosity of the GC via 
the following argument. In our models we are computing 
the scattered flux in response to an Eddington flare of Sgr 
A* some time ago, which enables us to estimate the ra- 
tio i^sc/-Pobs- Consistency with observations demand the 
scattered flux to be smaller than the actually observed 
one: i^sc/^obs < 1- Since the computed scattered flux just 
scales linearly with the assumed input luminosity we can 
find the luminosity of the GC for which this inequality is 
fulfilled: 

-^GC ^ -^-^Edd- (22) 

For twenty clouds selected from our sample we have listed 
the scattered flux, the observed flux and the derived limit 
upon the past luminosity of Sgr A* in columns 10-12 of 
Table 1. 

If the duration of a flare increases beyond td then the 

flux will stay at a constant level and the curves will extend 
to times t — to S ± A/2 as is illustrated in Fig. 4 by the 
responses to the flare lasting 300 years. An even further 
increase of the flare duration will have the result that at 
one point several clouds along one line of sight can be 
illuminated by one flare and their scattered fluxes begin 
to overlap, producing a larger total flux. This case will be 
discussed further below. 

As mentioned in the beginning of this section observa- 
tions show that the above assumed spherically symmetric 
density profiles are too simplified. The curves in Fig. 4 for 
a "real" GMC will probably look quite differently. They 
will be noisier, maybe show several peaks because of dif- 
ferent density maxima and will not display this symmetry. 
But, as we have shown, the choice of the density distribu- 
tion mainly determines the shape of the lightcurves, the 
scattered flux as a function of time, but not the level of 
the maximum scattered flux, which depends mainly on the 
total mass and the position of the cloud once it is fully cov- 
ered by a flare. We are therefore confident as long as the 
CO observations provide us with a good estimate of the 
total mass of a GMC that our computations, especially of 
the maximum scattered flux and the derived limit, are cor- 
rect, albeit the mentioned uncertainties about the density 
distribution inside a GMC. 

3.2. Response of one field of view 

The ASCA field of view has a diameter of about 40' and 
thus a size of 0.35 deg^. This diameter corresponds to a lin- 
ear scale of about 100 pc at the distance of the GC. Fig. 2 
shows that several clouds are located at similar Galactic 
longitudes. Their angular separations are often smaller 
than the diameter of the ASCA field of view. In most fields 
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Fig. 5. The scattered flux observable at a time t produced 
by five GMCs located in one ASCA field of view point- 
ing at a Galactic longitude of I = 27.51° in response to 
an Eddington flare of Sgr A* of 100 (dashed) and 1000 
(solid) years duration occuring around a time to- The ob- 
servational limit for the scattered flux produced in this 
field of view is marked by the dotted line. Note how the 
scattered fluxes of two GMCs begin to overlap for the 1000 
year flare because both arc illuminated. 



of view one finds therefore more than one cloud. Since we 
have a total of 96 GMCs distributed over 49 ASCA fields 
of view there are about two clouds per field of view on av- 
erage. For one of these fields of view (/ — 27.51°), in which 
five clouds are located, the maximum number of clouds in 
one field of view, we plot the scattered flux as a function 
of time in Fig. 5.^ 

Due to the different combinations of galactocentric and 
heliocentric distances of these five clouds they have differ- 
ent time delays, i. e. they respond to activity of the GC at 
different times. For every cloud we get a response similar 
to Fig. 4 at the time delay corresponding to each cloud. 

Fig. 5 shows that all five clouds already have reached 
their maximum fiux for the 100 year fiare. For this du- 
ration all the clouds are already covered and their satu- 
ration time is therefore smaller than 100 years. For the 
1000 year fiare the scattered fluxes from two clouds start 
to overlap. This happens when the duration of a flare is 
longer than the time needed for the ellipsoid to cross the 
distance between these two clouds. Despite the fact that 
Fig. 5 restricts the luminosity of Sgr A* much more than 



^ The plots presented in this and the following subsection 
have to cover a period of time of 40 000 years and also resolve 
flares on the timescale of years. In order to be more illustra- 
tive we are plotting the total scattered flux as a function of a 
different time coordinate defined by: 



t 



log 



We have chosen ta = 3000 years. This choice ensures a linear 
behaviour for times small and a logarithmic one for times large 
compared to ts- 



Fig. 4 it is obvious that there are large time gaps during 
the last 40 000 years where these 5 clouds are insensitive 
to the luminosity of Sgr A* . 

Spatial overlapping of the clouds in the fleld of view 
and absorption should not be a problem because the dense 
cores with large column depths lie in different directions. 
Certainly there should be further clouds of smaller size, 
below the detection threshold of the CO survey, present 
but these do provide only little scattered flux because of 
their much smaller masses. Since we have GMCs not only 
in one field of view one can combine the limits from all 
fields of view, which is done below. 



3.3. Response of all fields of view containing at least 
one GMC 

In total there are 49 ASCA fields of view covering at least 
one cloud out of the SRBY87 sample. The scattered fiux 
produced in these fields of view will be similar to the one 
shown in Fig. 5 with a peak for every cloud roughly pro- 
portional to the mass of the cloud. The peaks are centered 
around the time delay for the cloud and their temporal ex- 
tent is of the order of A/2, if the duration of the flare is 
longer than the cloud crossing time. 

To combine the information from all 49 fields of view 
we proceeded in the following way. We computed the ex- 
pected scattered flux for each field of view for a temporal 
grid covering the last 40 000 years with a sufficient time 
resolution to resolve the response of individual clouds. 
From these 49 lightcurves we select at each temporal grid 
point the largest scattered flux produced in response to the 
assumed input. As the observational limit at that time we 
choose the limit for that ASCA field of view, in which the 
largest scattered flux is produced. In Figs. 6-9 the result 
of this procedure is shown, the temporal response of all 
fields of view to an Eddington flare of Sgr A* of 1, 10, 100 
and 1000 years duration. 

As one progresses from short to long flare durations 
more and more time is covered by the GMCs. Increasing 
the flare duration by an order of magnitude results in an 
increase of the scattered flux by approximately the same 
factor (this factor depends upon the density proflle) as 
long as the flare duration is smaller than the "saturation 
time" for a given cloud. 

The first spike is produced by a GMC having a time de- 
lay of about 335 years (cloud 014 in Table 1). Remarkable 
about this cloud is that it attains its maximum scattered 
flux for a fiare duration of only about 12 years. With its 
radius of 5.8 pc one can estimate that the ellipsoid sweeps 
with a "velocity" of about 3.2 c across the cloud. This very 
large velocity is the result of the position of this GMC. It 
lies nearly half way between the GC and the observer at 
an longitude of only 10°. In Sec. 3. 1. we showed that 
clouds near the line between the GC and the observer can 
have very small cloud crossing times, because the motion 
of the ellipsoid is highly superluminal. 
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Fig. 6. The maximum of tlie scattered flux observable at 

a time t in all the fields of view produced by all 96 GMCs 
from our sample in response to an Eddington flare of Sgr 
A* of 1 year duration occuring around a time to- The 
observational limits for the scattered flux are marked by 
a dotted line. 
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Fig. 7. The maximum of the scattered flux observable at 

a time t in all the fields of view produced by all 96 GMCs 
from our sample in response to an Eddington flare of Sgr 
A* of 10 years duration occuring around a time to- The 
observational limits for the scattered flux are marked by 
a dotted line. 
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Fig. 8. The maximum of the scattered flux observable at 

a time t in all the fields of view produced by all 96 GMCs 
from our sample in response to an Eddington flare of Sgr 
A* of 100 years duration occuring around a time to- The 
observational limits for the scattered flux are marked by 
a dotted line. 
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Fig. 9. The maximum of the scattered flux observable at 

a time t in all the fields of view produced by all 96 GMCs 
from our sample in response to an Eddington flare of Sgr 
A* of 1000 years duration occuring around a time to- The 
observational limits for the scattered flux are marked by 
a dotted line. 



The largest scattered flux is produced by a GMC with 
a time delay of 669 years (cloud 003 in Table 1). To 
make the maximum scattered flux produced by this cloud 
(~ 7 X 10~* ergcm"^ s^^) consistent with the observa- 
tional limits (~ 2 X 10~^^ ergcm"^ s~^) the burst must 
have been ^ 3 x 10~^ times below the Eddington lumi- 
nosity, corresponding to a luminosity of Sgr A* of at most 
^ 8 X 10''°ergs~^. For the smaller peaks one derives less 
tight limits (see also Table 1). 

For short flare durations the flux from clouds with 
smaller time delays is larger than the one from clouds 
with large time delays. Increasing the flare duration these 



clouds attain comparable fluxes. This is mainly a conse- 
quence of the fact that the position of these clouds rela- 
tive to the GC and the observer, determines not only their 
time delay but also the time it takes the ellipsoid to cross 
the cloud. As we have seen the velocity of the ellipsoid is 
smallest for clouds "behind" the GC, which are also the 
clouds with large time delays. 

Although the times when there is no observable scat- 
tered flux shrink with increasing flare duration there are 
still some time gaps as shown in Fig. 9 even for an assumed 
flare duration of 1000 years. The two most pronounced 
time gaps lie around about 8000 and 14 000 years. It is ob- 
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Fig. 10. The fraction of time the scattered flux is larger 

than the observed upper limits during the last 40 000 
years as a function of the duration of the flare for 
different luminosities of Sgr A*: 10^" erg s~^ (triangles), 
10^^ erg s~^ (crosses) and 3 x 10^'* ergs" ^(squares). This 
quantity is an approximate estimate of the probability to 
detect a flare of duration A occuring during the last 40 000 
years. 



vious from the trend in Figs. 6-9 that by further increasing 
the assumed flare duration at one point the whole period 
of the last 40 000 years will be covered. This happens for 
flare durations longer than about 3000 years. This means 
that we can exclude a flare lasting longer than 3000 years 
of a luminosity larger than a few 10^^ ergs~^ to have hap- 
pened during the last 40 000 years. A similar conclusion 
can be drawn from Fig. 10. Here we tried to address the 
question of what is the chance that there was a strong 
flare of duration A during the last 40 000 years unnoticed 
by us? 

Fig. 10 shows for the clouds from the SRBY87 sam- 
ple the fraction of time the scattered flux was larger than 
the observational limits during the last 40 000 years for 
different flare durations. This quantity should be an ap- 
proximate estimate of the probability to detect a flare as a 
function of flare duration A. The curves should be contin- 
uous but we computed them only for several durations of 
the flare. It is obvious that this simple estimate does not 
take into account the level of the scattered flux above the 
observational limits. From this picture we conclude that 
an Eddington level flare of duration longer than about 
300 years occurring any time during the last 40 000 years 
is unlikely and one of duration longer than 3000 years can 
be excluded. Since as mentioned above we only were able 
to use the data from one Galactic quadrant the values of 
Fig. 10 actually should be regarded as lower limits to the 
probability. 



4. Limits on the X-ray luminosity of Sgr A* from 
the interstellar HI distribution 

Another major component of the ISM of the Milky Way 
is neutral, atomic hydrogen. Its total mass is about three 
times larger than that of molecular hydrogen. This is due 
to the larger extent of the HI disk. 

As shown in the previous section our sample of GMCs 
is insensitive to the activity of the GC during the last 
40 000 years for a few time windows, in particular two 
periods of 2000 and 4000 years duration of about 8000 
and 14 000 years ago for flares of about 1000 years dura- 
tion. For shorter flare durations these temporal gaps be- 
come larger and more numerous. They are as mentioned 
in Sec. 1 a result of the concentration of GMCs towards 
the molecular ring, their finite number in the Galactic 
disk, their finite size of about 10 pc and the probably the 
fact that we are missing the necessary data about GMCs 
located in the fourth Galactic quadrant. In this respect 
the more continuous HI distribution is of greater value. It 
provides less tight limits but extends further back into the 
past without large gaps in the temporal coverage. 

Interstellar neutral hydrogen is concentrated in clouds, 
which are much less massive (< 10^ Mq) than GMCs. 
Their total number is thus much larger. In the model of 
the Galactic ISM of McKee & Ostriker (1977) the mean 
separation between individual HI clouds A is of the order 
of 88 pc. Applied to our problem there exist two obvi- 
ous limits for the response of these clouds distributed in 
the Galactic disk to a flare of Sgr A* depending on its 
duration: 

a) If the fiare duration A is longer than the time for the el- 
lipsoid to travel the mean distance between HI clouds, 
A/c « 290 years, the flare scatters on a quasi continu- 
ous background medium. Because the gas is optically 
thin all the mass illuminated by the flare, the mass 
enclosed between the two ellipsoids, contributes to the 
scattered flux. 

b) For a flare duration shorter than the time it takes the 
ellipsoid to travel across the mean cloud separation 
we are dealing with individual clouds and the problem 
becomes more complex, in which each HI cloud gives 
some contribution to the total scattered flux. 

The HI distribution is rather thin with a scale height 
of about 100 pc (e.g. Dickey & Lockman 1990) for galacto- 
centric distances smaller than Rq. Beyond the solar circle 
the HI scale height seems to increase due to a possible 
warp of the Galactic disk. For distances larger than ~ Ro 
most of the disk is practically inside the ASCA field of 
view of 40'. For our purposes it is therefore valid to ne- 
glect the distribution of HI with height above the disk and 
to assume the HI to be concentrated in the Galactic disk. 
Since the HI mass surface density is approximately con- 
stant out to a radius of 16 kpc, we assume the following 
profile: 



^m{R) 



4 Mq pc-2 




(23) 
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Fig. 11. The scattered flux observable at a time t pro- 
duced by the Galactic HI distributed along one of the 
ASCA fields of view in response to Sgr A* switching off 
from Eddington luminosity at a time to. The different 
lines mark fields of view pointing at different Galactic lon- 
gitudes: ±10° (solid), ±30° (dashed), ±60° (dot-dashed) 
and ±90° (long dashed). The observational limit for the 
flux in the ASCA field of view pointing at I — 10° is 
marked by the dotted line. 

with i?max ~ 16 kpc. In this section we want to compute 
how the HI distributed in our Galaxy according to this 
profile responds to different types of activity of Sgr A*. 
The forms of activity we have in mind are: 

— Switch off: Sgr A* was active at a level LEdd up to a 
time to and then switched off. ^ 

— Flare: A flare lasting for a duration A at a level LEdd 
ended at a time to- 

For the Milky Way the switch off scenario is equivalent to 
a "flare" with a duration longer than about 110 000 years 
because the whole disk out to 16 kpc is illuminated in this 
case. 

For the mass surface density profile given by Eq. 23 the 
mass of hydrogen distributed along the ASCA beam point- 
ing at a longitude of 10° is of the order of 2 x 10'' Mq , which 
is about an order of magnitude larger than the most mas- 
sive GMCs in our sample. It is therefore of no surprise that 
if all of this gas is illuminated the scattered flux will be 
substantial. The maximal optical depth for Thomson scat- 
tering along an ASCA beam is of the order of 5 x 10~^. The 
main photoelectric absorption by heavy elements should 
be caused by GMCs along the way. One should therefore 
look for an ASCA field not "contaminated" by a massive 
GMC. 

In Fig. 11 we plot the scattered flux produced at a time 
t by the HI distributed along ASCA beams pointing at four 
different Galactic longitudes in response to a switch off of 
Sgr A* at a time to- The further away one looks from the 

* Note that in this and the following section to marks the 
end of the assumed activity. 




2x10* 4x10* 6x10* 8x10* 1x10^ 

t-t„/years 

Fig. 12. The scattered flux observable at a time t pro- 
duced by the Galactic HI distributed along one of the 
ASCA fields of view pointing at a Galactic longitude 
I = 10° in response to Eddington flares of Sgr A* of dif- 
ferent durations ending at a time switch off (solid), 
10 000 years (long dashed), 3000 years (dot-dashed) and 
1000 years (dashed). The observational limit for the flux 
in the ASCA field of view pointing at 10° is marked by 
the dotted line. 

GC the less mass will be located inside the field of view. 
This is the main reason for the differences of the scattered 
fluxes for different longitudes at times before the switch 
off. The cutoff of the curves at late times is due to the 
ellipsoid leaving the disk. For larger Galactic longitudes 
where the line of sight through the disk is smaller this 
occurs at earlier times, e. g. for 90° this happens after 
about 70 000 years whereas for 10° after 100 000 years. 

Fig. 12 shows the same quantities plotted for one 
Galactic longitude but for different durations of the flare. 
It is obvious that a switch off of Sgr A* from Eddington 
level less than about 100 000 years ago is inconsistent 
with the observational limits. For a luminosity typical of 
AGN hosted by a spiral galaxy of 10^'^ ergs~^ a switch off 
could not have happened less than about 80 000 years ago. 
Remarkable about this picture is that even a flare of a du- 
ration of 1000 years happening 100 000 years ago would 
produce a scattered flux larger than currently observed 
upper limits. A similar conclusion can be drawn from the 
following two figures. 

Fig. 13 shows the dependence of the scattered flux in 
an ASCA field of view at different times after the end of 
the activity in the longitude range —60° < I < 60°. Due 
to our assumed axially symmetric mass surface density 
profile the curves are obviously symmetric in the Galactic 
plane about the direction towards the GC. At small longi- 
tudes the scattered flux is large because one looks through 
a lot of gas and also the distance of the gas to Sgr A* is 
smaller than for larger longitudes. The scattered fluxes 
are rather large, because for the switch off case the whole 
disk is scattering at early times after the demise of the 
nuclear activity. The cutoff observable for the scattered 
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Fig. 13. The scattered flux produced by the Galactic HI 

distributed along one ASCA field of view in response to 
Sgr A* switching ofi^ from Eddington luminosity at a time 
to as a function of Galactic longitude I. The diff'erent lines 
mark the scattered flux observable at different times after 
the switch off of Sgr A* : 1000, 3000, 10 000, 30 000, 50 000, 
80000 and 100 000 years (from top to bottom). The ob- 
servational ASCA limits are marked by crosses. 

flux 100 000 years after the switch off at Galactic longi- 
tudes |/| > 35° is again a result of the ellipsoid growing 
beyond the outer edge of the HI disk and therefore leaving 
no illuminated matter along that direction. 

For a shorter flare the flux will be obviously less, be- 
cause only the illuminated parts of the disk will scatter the 
incident radiation. In Fig. 14 we plot the expected flux for 
a flare of 300 years duration. As argued above flares of 
this or a longer duration should scatter on a "continuous" 
HI background. 

The dependence of the flux scattered by the Galactic 
HI upon Galactic longitude in response to a flare with a 
suitable choice of duration, time elapsed since the fading 
of the flare and luminosity is similar to the X-ray emission 
of the Galactic ridge as observed by ASCA. This scenario 
might therefore be a possible explanation for the hard 
component of the X-ray emission of the Galactic ridge. 
Under the assumption that the flare of Sgr A* had a hard 
enough spectrum the scattered emission might contribute 
to the hard component of the emission of the Galactic 
ridge detected above 10 keV by RXTE (Valinia & Marshall 
1998) and OSSE up to energies of 600 keV (Skibo et 
al. 1997). There is one crucial test for this model. If we re- 
ally observe scattered radiation from Sgr A* then a narrow 
iron Ka-line at 6.4 keV must be present with an equiva- 
lent width of the order of 1 keV (Vainshtein & Sunyaev 
1990). 

The detection of emission in the 6.7 keV line from 
the Galactic ridge by Tenma was reported by Koyama 
et al. (1986). This fact shows that the upper limits from 
ASCA we are using are coming from the detection of the 
non scattered component. The detection of the neutral 
iron Ko,-line intensity in the same direction might give us 
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Fig. 14. The scattered flux produced by the Galactic HI 

distributed along one ASCA field of view in response to 
an Eddington flare of Sgr A* of 300 years duration ending 
at a time to as a function of Galactic longitude I. The dif- 
ferent lines mark the scattered flux observable at different 
times after the end of the flare: 1000, 3000, 10 000, 30 000, 
50 000, 80 000 and 100 000 years (from top to bottom). 
The observational ASCA limits are marked by crosses. It 
is possible for the flare to reach large longitudes because 
it moves with highly superluminal motion at early times. 

much stronger limits than presented in this paper. XMM- 
Newton is certainly able to separate the 6.4 and 6.7 keV 
lines much better than Tenma. Therefore we will get this 
information in the coming years and will be able to im- 
prove on the results of this paper. 

Obviously spiral arms or other large scale structures, 
a warp of the Galactic disk, will change the simple pic- 
ture described above. It will lead to increased scatter- 
ing for some directions, e. g. for large column densities 
along spiral arms, and thereby provide additional infor- 
mation about these deviations from a constant disk. If 
for example the density contrast of HI inside and outside 
of spiral arms is substantial one would get an additional 
time dependence as the ellipsoid scans across spiral arms 
and the regions between them. Nevertheless in this paper 
we decided to assume the simplifled distribution given by 
Eq. 23. 

5. The Milky Way at extragalactic distances 

It is well known that other spiral galaxies contain an ISM 

similar to the one of the Milky Way. CO and 21 cm ob- 
servations have revealed the spatial distribution of H2 and 
HI in many nearby spiral and irregular galaxies. M31 for 
example, containing a supermassive black hole of about 
4.5 X 10^ Mq in its center (see e. g. Kormendy & Gebhardt 
2001), has a HI mass similar to the one of the Milky Way, 
which is very concentrated in a ring at about 12 kpc from 
the nucleus (Sofue & Kato 1981). Although its H2 mass is 
less than the one of the Milky Way it is also located in a 
ring at about the same radius (Nieten 2001; Wielebinski 
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Fig. 15. The surface brightness of the scattered emission 
in the plane of the sky (specified by x and y or p and tp) 
produced by a fully illuminated gas disk of constant mass 
surface density 10 000 years after the central source has 
turned off. The inclination angle of the disk is i = 60°. 
The physical radius of the disk is 16 kpc. The center of 
the circle (marked by a cross) which marks the boundary 
between the illuminated and not illuminated parts of the 
disk moves with a velocity of 1.73c in the direction given 
by ?/i = TT. Its radius has a size of about 6 kpc and grows 
with time at a rate 2c. 



priv. comm.). Since there have been and will be Chandra 
(Garcia et al. 2000) and XMM-Newton (Shirey et al. 2001) 
observations of this object, M31 might be a good target to 
apply to the methods and results presented in this section. 
Many other nearby spiral galaxies have been observed by 
Chandra already, e. g. MlOO (Kaaret 2001), MlOl (Pence 
et al. 2001), NGC 1068 (Young et al. 2001) and NGC 4151 
(Yang et al. 2001) to name a few. 

Many of these spiral galaxies also seem to harbour a 
supermassive black hole at their centers (Kormendy & 
Gebhardt 2001) with some of them demonstrating AGN 
activity. Spiral galaxies are believed to be the host galax- 
ies of certain classes of AGN, in particular Seyfert I and 
II galaxies. The luminosities of their nuclei are typically 
of the order of 10^^ — lO^'^ergs"^. It is therefore an in- 
teresting question to ask how interstellar gas responds to 
activity of the AGN and what observational signs it should 
produce when the AGN switches off. We want to address 
these questions using the methods presented in the previ- 
ous sections. 

To demonstrate solely the effects of scattering by a gas 
disk we decided not to model a specific galaxy with an ob- 
served density distribution of its neutral gas but rather to 
show how the scattered emission of the Milky Way would 
look like at extragalactic distances. We therefore assume 
the same simple mass surface density profile as in the pre- 
vious section (Eq. 23). 



5.1. Scattered X-ray surface brightness 

If a galaxy is close enough to be resolved by the Chandra 
X-ray telescope then it is possible to map its scattered 
surface brightness as a function of position on the sky. 
The basic formula to compute this quantity along a given 



line of sight is Eq. 11. In this section we are using the 
paraboloid approximation (see Sunyaev & Churazov 1998 
or Sec. 1) to compute the integration boundaries Si„in(i) 
and Smax(i)- The paraboloid approximation is valid in this 
case given the large distances of extragalactic objects. 

In Fig. 15 we plot the surface brightness of the 
scattered emission produced by a gas disk with a cen- 
tral, isotropic radiation source of luminosity Lagn — 
lO''^ erg s~^ that turned off 10 000 years ago. The disk has 
an inclination oi i = 60°, with i defined as the angle be- 
tween the line of sight and the normal direction of the disk. 
The figure shows that a part of the disk is not bright, be- 
cause it is not illuminated anymore^. The boundary of this 
region, which is the projection of the intersection of the 
paraboloid at a time t — to with the disk onto the plane 
of the sky, is a circle. As time progresses this boundary 
moves radially outwards from the nucleus with an appar- 
ent velocity 



ccosi 



-I- COS ■0 sin i 



(24) 



This velocity depends upon the position angle ipj which is 
defined in Fig. 15. The maximum velocity of ccos«/(l — 
sini) is given for ip = tt and the minimum velocity of 
ccosi/(l + sini) for ip = Q. This "motion" has the effect 
that the center of the circle (marked by a cross in Fig. 15) 
moves in the direction -0 = while it is growing in size. 
The velocity of its center is given by Wq = (p(''i')~p(0))/2 = 
c tani, which is superluminal for inclination angles larger 
than 45°. For an inclination of i = 60° this velocity is 
Vo w 1.73c. The radius of the circle grows with time as 
ct/ cos i. 

The motion of this circle and the boundary between 
the illuminated and not illuminated parts of the disk is 
an example of apparent superluminal motion produced 
by scattering on suitable shaped screens. This subject 
has been discussed in the literature by e. g. Blandford 
et al. (1977). 

For a face on disk (i — 0°) with a constant mass sur- 
face density it is easy to see from Eq. 11 that the surface 
brightness drops with distance from the nucleus as cx p^^ . 
This is still valid in the case of an inclined disk. In the limit 
of a uniform disk with a constant mass surface density re- 
stricted by r = Tmax one can derive an analytic solution 
for the surface brightness of the illuminated parts of a 
disk with inclination i. It is necessary to express the fac- 
tor (1-1- cos^ dsc)/^^ appearing in Eq. 11, where O^c is the 
scattering angle and r the physical distance from the nu- 
cleus, as a function of the projected radius p, the position 
angle ip and the inclination angle i which yields 



^HI 



n 3(Tt 1 + 2 cos^ -0 tan^ i 
cos i Anp'^ 167r (1 -|- cos^ ip tan^ z)^ 



(25) 



As visible from Fig. 15, beyond a certain position an- 
gle the scattered surface brightness cuts off because these 



^ Note that the "negative" of Fig. 15 corresponds to the case 
10 000 years after a switch on of the AGN. 
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Fig. 16. The scattered X-ray surface brightness as a func- 
tion of projected distance from the nucleus along a fixed 
position angle oiip = tt/2. The solid line marks the surface 
brightness for the fully illuminated disk. The other lines 
mark the switch off (Fig. 15) and the 1000 year flare case 
(Fig. 19) 5000 (dot-dashed), 10 000 (dashed) and 20 000 
(dotted) years after the end of the nuclear activity. The 
switch off case corresponds to the curves which show no 
brightness below a certain distance and then join on to 
the solution for the fully illuminated disk. The increase 
in surface brightness between 4 and 6 kpc is due to the 
molecular ring which is located at these radii. The 
dependence described by Eq. 25 is visible. 

parts of the disk are not illuminated anymore by the cen- 
tral source. This angle is given by: 

V'max = arccos : . (26) 

\2ctpiBxii J 

Besides this more qualitative picture we present in 
Figs. 16, 17 and 18 the surface brightness along several 
cuts of the partly illuminated disk shown in Fig. 15. These 
figures confirm the validity of the above formulae, Eqs. 25 
and 26. 

In Fig. 16 we plot the scattered surface brightness as 
a function of projected radius p for a fixed position an- 
gle of -0 = tt/1. The curve for the fully illuminated disk 
shows the p^'^ dependence of the surface brightness given 
by Eq. 25. Note that in this computation we have included 
a molecular ring with a two times larger mass surface den- 
sity between about 4 and 6 kpc, which causes the bump 
between these radii. The surface brightness cuts off beyond 
the maximal disk radius of 16 kpc. Furthermore plotted 
in this figure are the surface brightness for the switch off 
(full picture given in Fig. 15) and the flare scenario (full 
picture given in Fig. 19) at different times after the end 
of the nuclear activity. The curves showing the surface 
brightness for the the flare scenario move with a velocity 
of about c to larger projected distances as given by Eq. 24 
for a position angle of V' = 7r/2. 

Figs. 17 and 18 show the scattered surface brightness 
at a fixed projected distance of p = 2 kpc from the nucleus 



Fig. 17. The scattered X-ray surface brightness as a func- 
tion of position angle at a fixed projected distance of 2 
kpc from the nucleus. The different lines mark different 
inclination angles of the disk for the fully illuminated disk 
and 5000 years after switch off: i ~ Q° (dotted), i = 30° 
(dashed), i = 60° (dash-dotted) and i = 80° (solid). The 
switch off leads to to disappearance of the illuminated re- 
gion at large position angles. 

as a function of position angle. For clarity no molecular 
ring was assumed in these two plots. 

In Fig. 17 we plot the surface brightness of the scat- 
tered emission for different inclinations of the disk for the 
fully illuminated and the switch off case. The brightness 
contrast for large inclinations of the disk mainly comes 
about because the same projected distance at position 
angles V = and V = 7r/2 corresponds to physical dis- 
tances from the nucleus differing by a factor of (cosi)"""^. 
Furthermore the average scattering angle for material at 
a position angle of V' = is ~ ''"/2 -|- i compared to 
Osc ~ 7r/2 for a position angle of ijj = tt/2. The angu- 
lar dependence of the curves presented in Fig. 17 is well 
described by Eq. 25. 

Fig. 18 shows the same quantities plotted for one disk 
with a large inclination oi i = 80° at different times after 
the switch off of the central source. With increasing time 
the not illuminated part of the disk "works" its way to- 
wards smaller position angles. The angle where the bright- 
ness cuts off is given by Eq. 26. 

In Fig. 19 we show how structure in the same disk 
as in the previous figures is probed by a fiare of the nu- 
cleus. This case might be relevant for spiral galaxies like 
Andromeda, where the neutral gas is very concentrated 
in a ring. We map the surface brightness of the scattered 
emission produced by a gas disk with a 2 kpc broad molec- 
ular ring at 5 kpc, with a mass surface density two times 
larger than the other parts of the disk, illuminated by a 
flare of the central source of 1000 years duration at differ- 
ent times after the end of the flare. The velocity of the light 
front behind (or above) the nucleus is much smaller than 
the superluminal one in front (or below) of the source. 
As time progresses the light front scans along the molec- 
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Fig. 18. The scattered X-ray surface brightness as a func- 
tion of position angle at a fixed projected distance of 2 Icpc 
from the nucleus. The inclination of the disk is « = 80°. 
We plot the surface brightness for the fully illuminated 
disk (solid line) and 1000, 3000, 10 000, 30 000 and 50 000 
years (dotted lines from right to left) after the switch off 
of the central radiation source. 

ular ring. GMCs will lighten up and and produce flares in 
the response. With the spatial resolution of Chandra one 
should be able to resolve these clouds in a galaxy up to a 
distance of about 10 Mpc. The characteristic time of these 
events is rather short. 

At a certain point the lower, faster light front reaches 
the edge of the gas disk. The "ring" will start to break up 
and a luminous arc travels upwards, slowly fading away 
as it reaches larger distances. 

Obviously a flare due to the tidal disruption of a star in 
the nucleus with a duration of the order of a few years will 
produce a feature in the scattered surface brightness dis- 
tribution of smaller width than the one shown in Fig. 19. 
If the tidal disruption rate is large enough, of the order 
of 10~*year~^, one might find several of these arcs per 
galaxy. 

5.2. Scattered X-ray luminosity 

For a source too distant to be resolved by X-ray telescopes 
the measurable quantity is the scattered flux from which 
one can derive a scattered luminosity. In this section we 
compute the scattered luminosity of a spiral galaxy with 
the same neutral gas distribution as our Galaxy in re- 
sponse to the above assumed forms of activity of its nu- 
cleus. 

Fig. 20 shows the scattered luminosity after the AGN 
has switched off as a function of time for different in- 
clinations of the disk. As parameters we choose as be- 
fore iagn = lO^^ergs"^ and a constant HI distribution 
Shi = 4M0pc^^ up to 16 kpc. 

The different luminosities for different inclinations at 
times t — to < 0, when the whole disk is illuminated, come 
about because of the angular dependence of the differ- 
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Fig. 19. The surface brightness of the scattered emission 

produced by a gas disk with a 2 kpc broad molecular ring 
at 5 kpc illuminated by a flare of the central source of 1000 
years duration having occured 5000, 10 000, 20 000, 40 000 
and 60 000 years ago. The edges of the gas disk at 16 kpc 
and the the molecular ring at 4 and 6 kpc are marked by 
solid ellipses. The inclination angle of the disk is i = 60°. 
The centers of the "rings" 5000 and 10 000 year after the 
end of the flare have been marked by a plus sign and a 
cross. Note how with progressing time the flare scans along 
the molecular ring. The highest surface brightness regions, 
besides the one close to the nucleus, are where the flare 
intersects with the molecular ring. 

ential scattering cross-section. For a face on disk i = 0° 
the scattering angle is obviously tt/2. Because the differ- 
ential cross-section for Thomson scattering is minimal for 
an angle of 7r/2 the face on disk has the smallest scattered 
luminosity. The much stronger differences at larger times 
are due to the super luminal motion of the light front. It 
leads to the strong effects in the low luminosity tail. Note 
the similarity of these curves with the ones of the scattered 
flux produced by the Galactic HI distribution in response 
to a switch off or a flare of Sgr A* (Figs. 11 and 12). 
The edge on case {i = 90°) is similar to the case of the 
observer located inside the Milky Way looking at small 
Galactic longitudes. The X-ray binary population of the 
Milky Way produces a luminosity of about 2 x 10^^ erg s"-'^ 
in the 2-10 keV range (Grimm et al. 2001). It will be there- 
fore difiicult to detect the scattered signal below this limit. 

For the face on disk [i — 0°) one can derive a simple 
analytic solution for the scattered luminosity as a func- 
tion of time. In this case the illuminated and not illumi- 
nated parts are concentric disks (switch off or switch on 
scenario) or annuli (flare scenario). The problem there- 
fore is axially symmetric and the projected distance p is 
equivalent to the physical distance r from the nucleus. 
Imagine the disk to made up of small mass elements 
dM = E(/9) d^ = S(/9)27rpd/0. If the mass element is il- 
luminated it contributes a fraction dL oc dM/p^ (Eq. 14) 
to the scattered luminosity. The total scattered luminos- 
ity at a time t is just the integral over all the illuminated 
mass elements 

L.M-f'iLo:rm±^ (27) 

J Ja{t) P 
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Fig. 20. The scattered X-ray luminosity observable at a 

time t produced by a spiral galaxy with the HI distribu- 
tion of the Milky Way in response to its central radiation 
source switching off at a time to- The assumed luminos- 
ity of the AGN is lO'^'^ergs-^ The different curves are 
for an inclination angle of the disk of 90° (solid), 60° (dot- 
dashed), 30° (long dashed) and 0° (dashed). The dotted 
line marks the Galactic luminosity produced by X-ray bi- 
naries. 

For a constant mass surface density S(p) = Shi one ob- 
tains the simple solution 

L,S) = L^:^^^y (28) 
with 

C = 2XI0»=„S.- (jjjJ^^) (j^) ^ 

For the face on disk the ellipsoids travel with the velocity 

of light radially outward and the the inner and outer ra- 
dius for the flare scenario can be be described as a{t) = ct 
and b{t) = c{t + A). For the switch off scenario one has 
b{t) = pinax and for the switch on case a{t) = 0. This sim- 
ple analytic formula helps to understand the behaviour of 
the curves presented in Fig. 21, which are the results of 
our numerical integration. For large times the flare solu- 
tions join with the switch off solution. This happens when 
only the outer parts of the disk of size corresponding to 
the flare duration are still scattering. 

5.3. Fluorescent iron Ka-line luminosity 

Scattering of X-rays produced by a powerful AGN is ac- 
companied by X-ray photo absorption by heavy elements 
and emission in the narrow fluorescent Kc-lines. The iron 
Ka-line at 6.4 keV is the brightest spectral emission fea- 
ture in the scattered radiation. Its equivalent width in the 
scattered X-ray flux is close to 1 keV for a normal cosmic 
abundance of iron (Vainshtein & Sunyaev 1980). The time 
dependence of the K„-line luminosity should repeat that 
of the scattered continuum. The main difference is that 
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l-t„/years 

Fig. 21. The scattered X-ray luminosity observable at a 

time t produced by a face on disk {i = 0°) with the HI 
distribution of the Milky Way in response to flares of the 
central radiation source of different durations ending at 
a time to: 1000 years (dashed), 3000 years (long dashed) 
and 10 000 years (dot-dashed). The solid line shows the 
switch off case. The assumed luminosity of the AGN is 
lO'''^ erg s~^. The dotted line marks the Galactic luminos- 
ity produced by X-ray binaries. 

scattering in the Ka-line is isotropic. The results for the 
K„-line fluorescence line should therefore differ by a factor 
of two or less and follow similar trends as presented in the 
previous figures. For example if we would plot the lumi- 
nosity in the K^-line in Fig. 20 the differences for times 
before the switch off would disappear and one would get 
the same luminosity for all disks independent of inclina- 
tion. Since the narrow Ka-line originating in Galactic X- 
ray binaries has a very low equivalent width, observations 
of the fluorescent, iron Ka-line are a more promising way 
to detect the scattered radiation produced by a faded flare 
of Sgr A*. 

6. Discussion 

In this paper we have presented upper limits on the hard 

X-ray luminosity of Sgr A* in the recent past. As the scat- 
tering material we have considered molecular and neutral 
atomic hydrogen present in the Galactic Disk. The CO 
and X-ray data we had access to, have allowed us to de- 
rive limits for the 4-10 keV luminosity of Sgr A* down to 
about 8 X 10^°ergs~^ at certain times. For other times 
the limits are of the order of 10''^ - 10''^ ergs . Now we 
want to briefly address a few points which demand some 
further discussion: 

a) Anisotropy of the X-ray emission of Sgr A* 

In simulating the activity of the GC we made the 
assumption that the emission of Sgr A* was isotropic. 
Studies of AGN show that their radiation field is often 
anisotropic due to the presence of a relativistic jet 
or as a result of shielding by a molecular torus. A 
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jet produced by Sgr A* would probably point away 
from the Galactic plane. This would make our derived 
limits weaker. 

b ) Chandra observations of the Galactic Ridge 
Ebisawa et al. (2001) have reported the results of a 
deep Chandra observation of a Galactic plane region 
(/ = 28.45°, 6 = -0.2°), which is devoid of known X- 
ray point sources. The total X-ray flux in the 2-10 keV 
range is determined to be ^ 1.1 x 10"^" ergs~^deg~^. 
About ten per cent of this flux is accounted for 
by point sources resolved by Chandra. These are 
extragalactic sources seen through the Galactic disk. 
90 % of the observed flux is due to diffuse emission at 
the level of the sensitivity and the angular resolution 
of Chandra. This diffuse flux is within a factor of 
two of the ones measured by ASCA in nearby fields, 
which are of the order of ~ 6 x 10~^^ ergs~^deg~^ 
(Sugizaki et al. 2001). Unfortunately none of the 
GMCs hsted in the SRBY87 sample lies in this field, 
but the HI column density in this direction is about 
2 X lO^^cm"^, similar to the one in the ASCA fields. 
Improvements might result from the spectral data 
obtained by XMM-Newton. If XMM-Newton will be 
able to detect the Ka-line of neutral iron in the same 
field this will open the possibility to separate the 
scattering contribution. 

c) Chandra observations oj the Orion Nebula 

The Trapezium region (/ = 209.01°, 6 = -19.38°) 
of the Orion Nebula has recently been observed by 
Chandra (Schulz et al. 2001). An upper limit for the 
diffuse emission in the energy range 0.1-10 keV of 
2 X 10^^ergs~"'^arcsec~^ is derived, which corresponds 
to a diffuse flux of ~ 1 x 10~^ ergs~^deg~^. The 
distance towards the Orion nebula is ~ 440 pc, 
corresponding to a time delay of ~ 2400 years. The 
galactocentric distance of Orion is i? ~ 8.8 kpc. 
The measured H2 column density in this direction is 
3.2 X 10^^ cm~^, which translates into an upper limit 
for the luminosity of Sgr A* of 3 x 10'*^ergs~^. So 
GMCs in our vicinity might give comparable limits 
to the ones derivable from GMCs located in the 
molecular ring, but their time delays will be only of 
the order of their distances. 

An obvious next step would be to apply this method to 
constr ain the luminositv of Sgr A* in the past to GMCs 



for XMM-Newton observations of the iron Ka-line in the 
direction towards massive gas complexes in the LMC. 
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